On the energetic origin of self-limiting trenches formed around Ge/Si quantum dots 
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At high growth temperatures, the misfit strain at the boundary of Ge quantum dots on Si(OOl) is 
reheved by formation of trenches around the base of the islands. The depth of the trenches has been 
observed to saturate at a level that depends on the base-width of the islands. Using finite element 
simulations, we show that the self-limiting nature of trench depth is due to a competition between 
the elastic relaxation energy gained by the formation of the trench and the surface energy cost for 
creating the trench. Our simulations predict a linear increase of the trench depth with the island 
radius, in quantitative agreement with the experimental observations of Drucker and coworkers Q. 



Heteroepitaxial growth of strained semiconductor thin 
films has been a subject of interest in the recent years, as 
it provides a versatile route for fabricating nanostructures 
that have potential applications as electronic, optoelec- 
tronic or memory devices. Lattice mismatch is one of the 
key elements that has been exploited to grow structures 
such as quantum dots and nanowires. Strain-driven self 
assembly in Ge/Si(001) has received particular attention, 
as 3D islands formed in this system can be integrated as 
devices that are compatible with well-developed Si tech- 
nology. It is well known that mismatch strain at low Ge 
coverages is relaxed by the formation of pyramid-shaped 
islands with {105} facets. As these islands grow in size, 
they transform into what are called domes with steeper 
sidewall angles. Even larger islands show strain relax- 
ation through the injection of misfit dislocations. In high 
temperature growth (T > GOO^C), strain relaxation is 
also achieved by trenches that form around the bases of 
both coherent and dislocated islands 

Experimental work of Chapparo and coworkers has 
shown that the trenches at the base of dome-shaped is- 
lands are self-limiting, in the sense that their depths sat- 
urate at a level, linearly proportional to the base- width of 
the island. While it has been suggested that this behav- 
ior is due to kinetic limitations [l| , quantitative estimates 
for the self-limiting depth and a reason for its linear de- 
pendence of the island size have not been provided. In 
this article, we consider the energetics and strain relax- 
ation of islands surrounded by trenches using finite ele- 
ment simulations and show that the self-limiting nature 
of trenches can be explained on the basis of a compe- 
tition between the elastic energy gained by the forma- 
tion of the trenches and the surface energy cost involved 
in creating the trenches. Our calculations that employ 
trenches whose dimensions and shapes are comparable 
to the ones observed in experiments, predict a linear de- 
pendence of the trench depth on the island base-width, 
in good quantitative agreement with experimental data. 

We first consider strain relaxation due to trench for- 
mation around a cone shaped island shown in Fig.^a). 
The cone radius and angle are denoted by R and 9, re- 
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FIG. 1: (a) Schematic of the axisymmetric model of a cone- 
shaped Ge/Si island of radius R and cone angle d. The depth 
of trench around the base of the island is denoted by h and the 
trench angles /3 and a are shown in the figure, (b) Contour 
plots of distribution of radial strain Err for an island with 
(right) and without trench (left). Here h = lOnm, R = 125, 
6=13 =25° and a=5°. 



spectively, while a trench of depth h is modeled using two 
inclined surfaces at angles (3 and a shown in Fig.lHa). In 
the case of trenches observed in experiments Jl, 2, 5], /3 is 
comparable to or bigger than the island sidewall angle 9, 
while the angle a is much smaller, usually a few degrees. 
A displacement based finite element method with linear 
shape functions is used to obtain the elastic fields in the 
island and the substrate. Strain in the structure can be 
expressed as 

where, e^j is the mismatch strain (present only in the is- 
land) and Uij is the additional strain due to elastic relax- 
ation. The non-zero components of the relaxation stress 
tensor (in cylindrical coordinates) are given by 
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where the elastic constants of the Si substrate and the Ge 
island are taken to be Cn = 165.7 GPa, C12 = 63.9 GPa 
and Gil = 128.5 GPa, C12 = 48.3 GPa, respectively. If 
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the strain fields are obtained from finite element calcu- 
lations, the free energy of the system, which is the sum 
of the energy gained through elastic relaxation and the 
increase in surface energy due to formation of the trench, 
can be written as 
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Here, the first integral is evaluated over the volume Vs of 
the system, which includes the island and the substrate 
with the trench, denotes the increase in surface area 
due to formation of the trench and 7 is the surface energy 
of the substrate, assumed to be independent of strain and 
orientation Q. 

The radial strain distribution in the island with and 
without the trench for (3 — 6 — 25", a = 5'^ and an 
equibiaxial mismatch strain of —4% is shown in Fig 1(b). 
In both the cases, the radial strain Srr changes sign across 
the island substrate interface with the substrate side un- 
der tension. The apex of the island is slightly over- 
relaxed The key difference between the case with 
the trench compared to the case without the trench re- 
lates to the strain distribution around the base of the 
island. The presence of trench leads to an effective re- 
laxation of the large near-surface compressive strain at 
point where the island makes contact with the substrate. 
At the same time, the trench also allows relaxation of 
tensile strain near the substrate surface. The azimuthal 
strain Ezz distribution has a maximum compressive value 
below the island boundary, so that the substrate near the 
boundary of the island tends to relax upwards and the 
material in the island tries to relax radially outwards. 

The elastic relaxation energy as a function of trench 
depth h is given in Fig.[21for islands of different radii, but 
fixed surface orientations 9, /3 and a. Here, we have plot- 
ted the first term in Eq. normalized by the total strain 
energy of the island prior to relaxation, Me^V, where M 
is the biaxial modulus = —4% is the mismatch 

strain and V — ttR^ tan 0/3 is the volume of the island. 
As the trench depth is increased, the elastic energy ini- 
tially decreases at a rapid rate, but eventually saturates 
at h/R ^ 0.3. Since the strain fields decay rapidly from 
the island substrate interface, the energy gained by ex- 
tending the trenches to relatively unstrained regions is 
small, leading to the observed saturation in strain relax- 
ation. We also find that normalized energies for islands 
of different radii fall on a single curve if the trench depth 
is scaled with the radius of the island, which shows that 
deeper trenches are required to achieve saturation of elas- 
tic energy for bigger islands. While elastic energy is re- 
laxed by increasing trench depth, deeper trenches also 
lead to an increase in the total surface energy of the sys- 
tem. The inset in Fig. 2 shows that the net change in 
surface energy is comparable to the energy gained via 
strain relaxation for an island with R = 125 nm. Fur- 
thermore, a competition between these two effects leads 
to a self- limiting trench size, ho, which corresponds to the 
minimum in the free energy F (refer to Eq. 2) as shown 



g -0.475 



Surface EnergyxlOO 
g -47.2 h 

1.00 ■ 



Free Energyxl 00 
■47.20 



o 

a 
DC 



-0.48 



-0.485 



-0.49 



-0.495 




R=100nm 

R=112.5nm 

R=125nm 

R=137.5nm 

R=150nm 



0.1 



0.2 



FIG. 2: Elastic relaxation energy of an island with a trench 
normalised by the total strain energy of the island prior to 
relaxation Me^V, plotted versus the ratio of trench depth h 
to island radius R. Note that the relaxation energy attains 
saturation at h/R ~ 0.3. The surface orientations (refer to 
Fig. 1) in this calculations are a= 5° and P=e^25°. The inset 
shows the competition between the elastic relaxation energy 
(blue) and the change in surface energy (green, normalised 
by M€qV) resulting in a minimum in the total free energy 
(black) at an optimum trench depth (ho) for an island with 
R = 125 nm and 7=1.5 J/m^ 

TABLE L Effect of variations in the surface orientations /3, a 
and 9 on the optimum trench depth {ho in nm) and the nor- 
malized relaxation energy {AF x 10"'^, where AF — F{h — 
0) — F{h — ho)) with respect to the ungrooved structure. 
When the variations of a particular parameter is considered, 
the default values of other parameters are, 9 — /3=25°, a=5° 
and 7=1.5 J/m^. 
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in the inset in Fig. 2. We point out that while the gain 
in elastic energy by formation of trenches is only 2% of 
the total strain energy of the structure without trench, 
this corresponds to 3.55xl0^fcBr at 600°C for the island 
in Fig. 2; the trenches should then be very stable against 
thermal fiuctuations. The net gain in free energy and the 
self-limiting trench depth depends on the surface orien- 
tations of the island and the trench and the size of the 
island. We consider these effects in the following para- 
graph. 

The data in the first column of Table ^ shows that 
the optimum trench depth ho increases by about Inm 
when the trench angle f3 is increased from 25° to 47" at 
fixed 9 and a for R — 125nm. We can understand this 
trend by noting that larger /? leads to more effective re- 
laxation of the compressive strain at the point where the 
island meets the substrate and the tensile strain in the 
substrate. While the free energy of the system can be 
lowered for larger /3, the opposite is true for the trench 
angle a (refer to the second column of Table P) . Lower 
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FIG. 3: Self-limiting trench depth versus island diameter for 
a conical Ge/Si island with 9=13— 25'\ q=5° and for the dome 
shaped island shown in the inset with H/R = 0.3 and /3=25°, 
0=5". The surface energy is 7=1.5 J/m^ in both the cases. 
Saturation depth observed in the experiments P| is shown by 
the filled circles. 

values of a allow better relaxation of the tensile strain 
in the substrate and involve a smaller change in surface 
energy. This result is in agreement with experimental 
observations Q, H, IE IE IB| that show steep (large P) and 
shallow (small a) surface orientations of the trench close 
to and away from the island, respectively. If the sidewall 
angle of the island is lowered by choosing j3 = 9 and by 
holding a constant, optimum trench depth /iq decreases 
in a monotonous manner (refer to the third column in 
Table. until it becomes unfavorable to form trenches 
for 6 < 15°. This result would indicate that the trenches 
would not be observed at around the hut-shaped {105}- 
faceted islands whose base- widths are « 200 nm, in agree- 
ment with experimental observations. Next, we consider 
the evolution of the trench depth with the radius of the 
island, with the aim of comparing the predictions of our 
calculations with the data reported in Ref. 1. 

Since the AFM and TEM images of the cross-sections 
of the domes 0,IEI3] indicate that the steeper face of the 



trench is at the same orientation as the sidewall of the 
dome cluster, we take (3 = = 25", which is very close 
to the angles made by the {113} and {102} facets with 
the substrate. The dependence of the optimum trench 
depth as a function of the island radius is given in Fig.[3| 
for the conical island in Fig. 1 and for a dome-shaped 
island shown in the inset in Fig. 3. The axially sym- 
metric sidewall of the latter shape is an arc of a circle 
with an aspect ratio (ratio of height H to radius R) of 
0.3. In both cases, we find nearly linear increase of the 
optimum trench depth with the island radius, which is 
in good quantitative agreement with experimental data 
from Ref. 1. It should be pointed out that our calcu- 
lations predict somewhat larger slopes for the curves, 
which could be attributed to a couple of modeling as- 
sumptions. In our calculations, we have taken that the 
composition of the islands is uniform (100% Ge), while 
there is experimental evidence that Si diffuses from 
the substrate into the island at temperatures considered 
in Ref. 1. In this case, since the elastic energy gain would 
be smaller than our estimates, the trench depths would 
also be somewhat smaller, bringing our results closer to 
the experimental data. Closer agreement could perhaps 
be achieved if strain and orientation dependence of sur- 
face energy is included in the analysis. 

In conclusion, we have shown that the self-limiting be- 
havior of the trenches around the bases of SiGe quantum 
dots can be explained on the basis of a competition be- 
tween elastic relaxation and the surface energy cost of 
creating the trench surface. The optimum trench depth 
has an approximate linear dependence with the island 
radius, due to the fact that more elastic energy can be 
relaxed by trench formation in the case of larger islands. 
Our results indicate that the saturation of the trench 
depths can be attributed to energetics rather than ki- 
netic limitations on mass transport. 

The research support of the National Science Founda- 
tion through grants CMS-0093714 and CMS-0210095 and 
the Brown University MRSEC program through grant 
DMR-0079964 is gratefully acknowledged. 



[1] S. A. Chaparro, Y. Zhang and J. Drucker, Appl. Phys. 

Lett. 76, 3534 (2000). 
[2] S. A. Chaparro, Y. Zhang, J. Drucker, D. Chan- 

drashekhar and D. J. Smith, J. Appl. Phys. 87, 2245 

(2000). 

[3] J. Drucker, Y. Zhang, S. A. Chaparro, D. Chan- 

drashekhar, M. R. Mccartney and D. J. Smith, Surf. Rev. 

Lett. 7, 527 (2000). 
[4] J. Drucker, IEEE J. Quantum Electron. 38, 975 (2002). 
[5] U. Denker, O. G. Schmidt, N. Y. Jin-PhiUipp and K. 

Eberl, Appl. Phys. Lett. 78, 3723 (2001). 
[6] U. Denker, M. W. DashieU, N. Y. Jin-PhiUipp and O. G. 

Schmidt, Mater. Sci. Eng. B. 89, 166 (2002). 
[7] T. I. Kamins, E. C. Carr, R. S. Williams and S. J. Rosner, 



J. Appl. Phys. 81, 211 (1997). 
[8] J. A. Floro, E. Chason, R. D. Twesten, R. Q. Hawang 

and L. B. Freund, Phys. Rev. Lett. 79, 3946 (1997). 
[9] Since data on the strain and orientation dependence of 
surface energies is not available from experiments or ac- 
curate quantum mechanical calculations, we use 7 = 1.5 
J/m^, which is typical for semiconductor surfaces. 

[10] A similar tensile strain at the top of islands has also been 
observed in calculations reported by H. T. Johnson and 
L. B. Freund, J. Appl. Phys. 81, 6081 (1997) and B. J. 
Spencer and J. Tersoff, Phys. Rev. B 63, 205424 (2001). 

[11] The biaxial modulus can be written in terms of the elastic 
constants as M = Cn -f C12 - 2C12/C11. 



